A wide family of singularity-free cosmological models 
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In this paper a family of non-singular cylindrical perfect fluid cosmologies is derived. The equation 
of state corresponds to a stiff fluid. The family depends on two independent functions under very 
simple conditions. A sufficient condition for geodesic completeness is provided. 
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I. INTRODUCTION 



After the discovery of the first regular perfect fluid cosmological model by Senovilla one of the major questions 
that that spacetime posed was to determine how generic this lack of singularities was. For instance, in the Ruiz- 
Senovilla ,2j barotropic family with equation of state p = 7/i, < 7 < 1, the regular models were those corresponding 
to a radiation fluid, 7 = 1/3. These solutions are the separatrix between those models with singular pressure and 
energy density, 7 > 1/3, and those who just had singular Weyl curvature, 7 < 1/3 Q. Therefore regular models 
would be a zero-measure set in this family. 

Our aim in this paper is to determine an infinite family of regular cosmological models that need not be so restrictive 
as the Ruiz- Senovilla family and therefore may indicate that regular models cannot be neglected in the set of solutions 
I of Einstein equations. 

5h . We shall focus on stiff perfect fluids, since they are simple enough to allow almost complete integration of Einstein 
^P" equations and thereby constitute an excellent arena for checking hypotheses. We shall show that under very simple 
^ , restrictions regular solutions appear. 

The second section of this paper will be devoted to deriving solutions of the stiff fluid Einstein equations in a 
convenient manner for our purposes. In the third section geodesic completeness of the solutions shall be imposed and 



the restrictions derived from this assumption shall be expressed as a sufficient condition. 



II. STIFF FLUID COSMOLOGIES 



We shall restrict to spacetimes endowed with an Abelian orthogonally transitive group of isometries, G2, acting on 
timelike surfaces, since this is the framework where regular cosmological models have so far appeared in. We shall 
further require that the Killing fields be mutually orthogonal. Adapting the coordinates to these fields we shall write 
them as {d z , d^}. Under these assumptions we can write the metric for the spacetime in a convenient coordinate chart 
0, {t,r,z,cl>}, 

ds 2 = e 2K (~dt 2 + dr 2 ) + e- 2U dz 2 + p 2 e 2U d0 2 , (1) 

which has been chosen as isotropic on the non-ignorable coordinates, t, r. This is a priori no restriction and may 
always be achieved since every 2-metric admits an isotropic parametrization. The metric is written in terms of three 
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functions, K, U and p that depend only on t and r. 

We may interpret the isometry group as cylindrical symmetry in the spacetime provided we have a regular axis 
where the norm of the angular Killing field vanishes. We shall come back to this issue later on. The range of the 
coordinates will be then the usual one for cylindrical symmetry, 

-oo < t. z < oo, < r < oo, < < 2ir. (2) 

The matter content of the spacetime will be a perfect fluid of energy density fi, pressure p and 4-velocity u. The 
energy-momentum tensor will be then, 

r"" = fivfu" +p{g"" + u»u v ), 0<M,f<3, u"u M = -l. (3) 

For a stiff fluid, (J, = p. We shall write down the Einstein equations, i?^„ — Rg^/2 — T^, in a comoving system 
of coordinates for the perfect fluid, that is, u = e~ K 3 t . After some simplications the equations read, 

Utt-U rr + -{U t pt-U r p r ) = 0, (4a) 
P 

Ptt - Prr = 0, (4b) 

K t p r + K r p t = p tr + U t p r + UrPt + 2pU t U r , (4c) 

K tPt + K rPr = Ptt % Prr + U tPt + U rPr + p{V 2 t + Ul) + ppe 2K , (4d) 
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U r p r - U t pt . TT 2 j T 2 2K 



K r + ^= 0, (4f) 



2p 

El 

p ' 2p 



K t + P - + £ = 0, (4g) 



where the last pair of equations are just Euler and continuity equations for the perfect fluid. 

Since every regular cosmological model in the literature has a p with spacelike gradient, we shall impose as an 
ansatz that grad p be orthogonal to the velocity of the fluid, u. 

Under this assumption, p will be a function of r only. But then l|4b(l requires that p be a linear function of r. After 
rescaling the coordinates we can take p — r and the whole system of equations becomes rather simple, 



U tt - U rr - — = 0, (5a) 
r 

K t = U t + 2rU t U r , (5b) 

K r = U r + r{U? + U?)+ pre 2K , (5c) 

K rr -K tt + ^ + U 2 - U 2 = pe 2K , (5d) 
r 

K r + 7T = 0, (5e) 



K t + E± = o. ( 5 f) 



2 P 
Pt 
2p 

The energy-momentum conservation equations can be integrated, 

p = ae- 2K , (6) 



with a = const. > and equation (|5dl) is a consequence of the others. We are left then with a two-dimensional 
reduced wave equation in polar coordinates without source term (|5a|l and a quadrature for K, 

K t = U t + 2rU t U r , (7a) 
K r = U r + r(U 2 + U 2 ) + oar, (7b) 

that can be integrated after providing a solution to the wave equation. The integrability condition for the quadrature 
(|7|) is namely the wave equation, so the whole problem reduces to solving it. 
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Following |5| , for instance, the solution to the Cauchy problem for the wave equation in the plane can be constructed 
from the 3D solution by ignoring the third variable, 

U(x,y,t) = ±- /% fdRR 9{* + R<x»*,V + R™*) , 1 8 fx I* fc + Rcnfiy + Rd**) 

for initial data U(x,y,0) = f(x,y), Ut(x,y,0) = g(x,y). Note that the 2D wave equation does not satisfy Huygens' 
principle and therefore the domain of dependence is a circle, not just a circumference. 

Our problem is even easier since we do not have dependence on the polar angle. Therefore we just have to impose 
circular symmetry on the initial data U{r, 0) = f(r), U t {r, 0) = g{r). The time coordinate may be removed from the 
integration limits by an appropriate scaling, R = tr, 

U(r, t) = ± fj dtf dr-^= \tg{v) + f(v) + tf(v) ^ + *J °° S * } , (9) 

where v = ^/r 2 + t 2 r 2 + 2rtr cos <fi, choosing the origin of the polar angle at the angle for (x,y). This expression is 
valid for all values of t. 

For instance, the non-singular spacetime in j(J is generated by 



f{x) = f * 2 , 9(x) = 0, (3 > 0. 

It is clear now that we are just integrating the functions /,j on a finite interval of r and that the integrals are 
well-defined provided g, /, /' are continuous. The singularity at t = 1 is harmless under such conditions. 

The solution does not share the class of differentiability of the initial data because of the derivative term. We need 
at least / £ C 3 ([0, oo)), g £ C 2 ([0, oo)) in order to have U,K £ C 2 ([0, oo)) and a well-defined Ricmann tensor. 

Surprisingly there is no need to impose cylindrical symmetry on the solution, since we have a regular axis at r = 
provided that /, g are regular there and we have already required it. According to Q, the axis is regular if 

< B ^A 8 nuiA) =ea(t7 _ K) , M =r 2 e 2 \ (10) 

r^0 4/\ 

therefore we need that K(Q, t) = U (0, t) for every value of t. But at the axis the equations I|7al7b|l that determine K 
are rather simple, 

K r = U r , K t = U t , (11) 

and therefore K(Q,t) — U(0,t) + const, and the condition of regularity at the axis is fulfilled by either taking the 
constant of integration equal to zero or conveniently rescaling the angular coordinate. 

Note that this requirement of regularity excludes a timelike gradient of p in the vicinity of the axis. 



III. GEODESIC COMPLETENESS 



The metric that is obtained after integrating the system (|7al7b|) has regular components in the whole spacetime, 
but this not suffices in order to have a non-singular spacetime. We shall consider that a spacetime is regular 
if it is causally geodesically complete, that is, if every causal geodesic may be extended to all values of its affine 
parametrization. 

This means analysing the geodesic equations for diagonal cylindrically symmetric spacetimes. This has been done 
in Those results can be summarized as follows: 

Theorem: A cylindrically symmetric diagonal metric in the form Q with C 2 metric functions /, g, p is future causally 
geodesically complete provided that along causal geodesies: 

1. For large values of t and increasing r, 

(a) K r + K t > 0, and either K r > or \K r \ < K r + K t . 

(b) {K + U) r + (K + U) t > 0, and cither {K + U) r > or \{K + U) r \ < (K + U) r + (K + U) t . 

(c) (K - U - \np) r + (K - U - Lnp)f > 0, and either (K — U — In p) r > or \{K - U - lnp) r \ < 
(K-U-ln p) r + (K-U-]xl p) t - 



4 



K(t,r)-U(t,r) 

2. For large values of t, constant b exist such that 2 K(t, r) ^ > — In |i| 

K(t,r) + U(t,r) + lnp(t,r) 

A similar result can be stated for past-pointing geodesies just reversing the sign of the time derivatives in condition 

1. 

Note that we have omitted the condition for non-radial geodesies with decreasing r in |8(, since according to (97) 
in that reference, non-radial causal geodesies would reach r — for finite t, that is, before t becomes singular, 
contradicting the fact that geodesies should be singular there. Therefore the axis cannot be reached by geodesies with 
non-zero angular momentum. 

All we have to do now in order to have a geodesically complete model is to check whether conditions 1 and 2 are 
satisfied. 

We first show that the conditions on the derivatives are always fulfilled for a stiff fluid model: 



1. (a) According to (|7al7b() we have, 

K t + K r = U t + U r + r(U t + U r j 2 + ar. 

Several possibilities are open: When Ut + U r is positive, K t + K r is positive. If Ut + U r is negative and 
\Ut + U r \ > 1, the quadratic term is larger and K t + K r is again positive for large values of r. Finally, if 
Ut + U r is negative and \U t + U r \ < 1, it is the pressure term ar which overcomes the negative term for 
large r. 

The same sort of reasoning is valid to conclude that the radial derivative, 

K r = U r + r(U? + U 2 r ) + ar. 

is positive for large values of r. 

(b) We can apply the same argument to K + U, 

(K + U) t + {K + U) r = 2(U t + U r ) + r(U t + U r f + ar, 

in order to show that these derivatives and (K + U) r are positive for large values of the radial coordinate. 

(c) Finally, the third condition is always satisfied, 

(K - U -In p) t + (K - U - lnp) r = r(U t + U r f - - + ar > 0, 

r 

(K-U- lnp) r = r(U? + U 2 ) -- + ar>0, 

r 

for increasing r. 

Past-pointing geodesies are treated analogously without any additional problem, since reversing the sign of time 
derivatives does not alter the positivity of the quadratic terms. 

Consequently only condition 2. yields a restriction. 

2. As we shall see, this condition amounts to study U at the axis for large values of the time coordinate: 

(a) This condition is trivial since (K — U)\ r= o = 0, 

(K-U)(t,r)= / dr'(K-U) r {t,r')= / dr' {r'(U 2 (t, r') + U 2 (t 7 r')) + ar} > 0. 
Jo Jo 

(b) The previous reasoning for ruling out singularities for decreasing radius leaves us with two possibilities: 
increasing radius and constant radius, r — 0. Since, according to 

K(r,t) = U(0,t) + [ dr'K r (r',t), 



and we have already checked that K r is positive for large t and increasing r, we just have to study the term 
U\ r= o, that is, we have reduced the problem to analysing the behaviour of K at the axis for large values 
oft. 
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(c) Similarly, 

(K + U) (r, t)+ In r = 2U(0,t) + In r+ [ dr (K + U) (/, i), 

Jo 

and we have already checked the positivity of K r + U r as in the previous condition. The logarithmic term 
does not mean a problem for increasing radius. Again we are left with controlling the behaviour of the U 
term. 

Summarizing our results so far, in order to have a causally geodesically complete spacetime we just have to require 
that U\ r= o does not decrease faster than a negative logarithm for large values of the absolute value of the time 
coordinate. The condition on the solution of the Cauchy problem for the wave equation at the axis becomes a bit 
simpler, 



K(0, t) = U(0, t) = J dr-^= {tg(\t\r) + f(\t\r) + \t\rf'(\t\r)} > -\ In \t\ + b, (12) 

since the dependence on the polar angle is lost. 

This bound can be attained, for example, by the initial data, 

\a\ - 1 a lnr 
f{r) = — 7, — ln r> 9(r) = , 

since for this choice of functions the solution to the Cauchy problem for the wave equation is, 

1 1 ( |(2|a| -l)ln|t|-ln2 sign at > 
U(0,t) = -(\a\ +asignt- l) ln]t| + -(|o| - asignt- l) ln2 = < 

2 2 [ | (2|a| — 1) In 2 — | In. jt| sign at <0 

The behaviour of the terms in (|12fl is rather different. The term, Uf , dependent on /, the initial value of U, is even 
in the time coordinate, as it is to be expected when the initial time derivative of U is zero. On the contrary, the term 
U g , dependent on g, is odd in t. 

This means that if Uf satisfies (|12|l for positive time, it is automatically satisfied for negative time. On the contrary, 
if U g satisfies (|12|l for positive time, it is only satisfied for negative time if it is also satisfied by U- g for positive time. 
Therefore three different possibilities are open depending on the value of 

hm m±im. (13 ) 

r^oo rg[r) 

• If lfH5|) is infinite, we need Uf{t) > — \ In \t\ + b for large values of t in order to have geodesic completeness. This 
means that f(r) + rf'(r) > — | lnr + k for large r. 

• If (fH5f is zero, we need |?7 9 (t)| < \ In \t\ + b for large values of t in order to have geodesic completeness. This 
means that \g(r)\ < lnr + k for large r. 

• If (|13|) is finite, then U / and U g are of the same order for large values of t or of —t and geodesic completeness 
will depend on the value of the limit. 

Note that values of the integrand close to r = do not influence the result for large values of t, since we may split 
the integral in two terms, 

fXc/t pi 

U(0,t) = / + / , 

J0 Jx a /t 

and the first one is bounded and negligible for large t. Therefore the main contribution to U comes from the second 
term, which must fulfill the required asymptotic behaviour. 
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IV. EXAMPLES 

A simple and wide family of functions that satisfy (|12|) can be written in terms of polynomials. Consider 



/(r)=J>r*, g(r) ( 14 ) 

i=Q i=0 

If n, m are even numbers, U can be analytically integrated in terms of polynomials. For our purposes we just require 
U at the axis, which can be integrated for a larger set of functions. Since U is linear in / and g, we may analyse the 
monomials independently. For f(r) = r™, g(r) = r m we obtain, 

n!! /7T\( 1 +(- 1 )" +1 )/ 2 m\\ /n\ (i+(-i) m+1 )/2 

•w-smd) u M-wmk) (15) 

These expressions are valid even for n,m = —1 taking (—1)!! := 1, although they may be not very practical. 
According to Ijl2|l we have two different possibilities for obtaining a singularity-free model: 

• If /, g are polynomials in r respectively of degree n,m and n > m + 1, we have a non-singular model if a n is 
positive. 

• If /, g are polynomials in r respectively of degree n,n — 1, Uf and U g at the axis are polynomials of degree n 
and we have a non-singular model if Uf dominates over U g . This happens if the leading term oiUf is greater 
than the one of U g , that is, 



2/ 2\ ( - 1} " Or, 



> 1. (16) 



>-l)!!y W |&n-l 

Using Stirling's formula for approximation of factorials, an easy and safe bound would be 

n+ an > |fc„-i|. (17) 

This family of non-singular cosmological models is large indeed, as it can be seen by restricting the range to 
a finite dimensional space of polynomial functions: if we consider the space of functions U for which U\ r =o 
is a polynomial of degree equal or lower than n, the subset of singularity-free models comprises an open set, 
according to fTTjl . 

This result can be generalized, since 

1 , 0F r(( P + 2)/2) 

rfT vT=^ = -r(b + 3)/2) (18) 

allows integration for every real value of the exponent p. Therefore, for f(r) — r p , g(r) = r q , we obtain, 
TTM /~ r((p + 2)/2) 0F r((g + 2)/2) 

which allow generalization of the geodesic completeness requirements which have been previously derived for 
polynomial functions to linear combinations of powers of r. Additionally one has just to care about the class of 
differentiability of U, which demands that p > 2, q > 2. 



V. CONCLUSIONS 



We have analysed a wide family of stiff perfect fluid cosmological models with cylindrical symmetry. The issue of 
causal geodesic completeness has been reduced to just the behaviour at the axis of the initial value problem for a 
sourceless 2D-wave equation, which is the only one left after simplifying Einstein equations. A sufficient condition for 
geodesic completeness is provided, which is very easy to check and to implement. The case of polynomial initial data 
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has been discussed and allows a fairly large set of non-singular cosmological models. We think that this set is wide 
enough to preclude considering non-singular models as isolated points in a space of cosmological models. 

The role of pressure in these models is obviously determinant, since stiff perfect fluids are a limit case for energy 
conditions, corresponding to a sound velocity equal to that of light. On the contrary, dust perfect fluids are always 
singular according to Raychaudhuri equation. Intermediate cases remain open for discussion, even though partial 
results have been obtained Q. 

It is interesting to notice that nonseparability of the models in these coordinates is fundamental for geodesic 
completeness. In |9( separable cosmologies were studied and none of them was found to be regular. 
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